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Introduction 


This  past  period  has  been  devoted  to  finalizing  the  lidar  system  design  and  ordering 
required  equipment.  A  great  deal  of  time  was  spent  on  studying  characteristics  and 
prices  of  intensified  CCD  systems.  Such  devices  fill  a  critical  need  since  they  are  to 
be  used  to  determine  angle  of  arrival  fluctuations  in  actual  images  of  the  lidar  returns 
from  different  elevations.  They  thus  provide  the  basis  for  measuring  the  elevation 
profile  of  C„2.  Such  a  device  is  also  the  most  expensive  component  of  our  planned 
lidar.  Consequently,  the  decision  as  to  which  model  to  purchase  from  which  company 
required  several  months  of  data  discussion  with  different  vendors.  In  the  end,  the 
Andor  ICCD  system  manufactured  by  Andor  Technology  was  chosen.  Finally,  all  the 
major  components  of  the  lidar  have  been  ordered,  and  most  have  already  been 
delivered.  The  intensified  CCD  should  arrive  within  two  months. 

In  addition,  a  new  lab  was  built  for  the  lidar.  The  room  is  located  on  the  roof  of  our 
department  building  and  is  characterized  by  two  windows  in  the  ceiling  -  one  above 
the  transmitter  and  one  above  the  receiver.  The  structural  changes  have  been 
completed. 

The  lidar  system  should  be  operational  within  a  month  or  two.  Even  though  the 
project  is  now  officially  completed,  it  is  our  intention  to  continue  issuing  progress 
reports  on  the  data  and  analyses  so  that  a  great  number  of  measurements  can  be  taken 
and  used  to  model  accurately  aerosol  size  distribution  and  turbulence  strength  vertical 
profiles  for  different  seasons  and  times  of  day,  and  their  effects  on  laser  beam 
propagation  at  various  altitudes. 

> 

These  should  provide  more  reliable  and  lasting  data  than  the  balloon  measurements  - 
suggested  in  the  original  proposal.  It  should  be  noted  that  the  balloons  and  Air  Force 
crews  to  make  the  originally.,  proposed  measurements  never  arrived.  The  lidar 
measurements  should  provide  more,  better,  and  more  reliable  data  than  the  originally 
proposed  balloon  measurements. 

The  lidar  system  and  expected  results  are  summarized  in  the  following  pages.  These 
expected  results  are  to  be  compared  with  the  actual  measurements.  We  will  not  be 
surprised  if  measurements  in  the  Middle  East  prove  to  be  somewhat  different  from 
those  elsewhere.  In  any  event,  as  shown  here  the  calculations  based  on  data  in  the 
literature  indicate  that  laser  beam  widening  at  elevations  from  2  to  20  km  caused 


by  aerosols  can  be  expected  to  be  very  significant,  quite  probably  even  more  so 
than  that  caused  by  turbulence.  Hence,  in  airborne  laser  weaponry  it  is  most 
prudent  to  consider  aerosol  effects  in  system  design  and  analyses  of  experiments. 


1.  Expected  lidar  returns. 


The  constituents  in  the  atmosphere  that  interact  with  the  transmitted  laser  pulse  to 
produce  the  backscattered  signal  include  the  gaseous  air  molecules  and  aerosol 
particles.  Molecules,  because  of  their  small  size,  scatter  light  rather  weakly,  producing 
Rayleigh  scattering  which  varies  inversely  with  the  fourth  power  of  wavelength  of  the 
incident  radiation  (actually  A,"40117).  Atmospheric  aerosol  particles  are  polydisperse 
and  span  a  wide  particle  size  range  (~0.01  to  10  micrometers)  which  includes  the  lidar 
wavelength  range.  Thus,  aerosol  particles  strongly  scatter/diffract  light  producing 
what  is  generally  approximated  as  Mie  scattering  after  the  generalized  spherical 
particle  scattering  theory.  At  a  given  lidar  wavelength,  molecular  scattering  varies 
directly  with  the  atmospheric  molecular  number  density,  but  aerosol  scattering  varies 
in  a  complicated  manner  depending  on  the  size  distribution  ,  shape,  and  refractive 
index  of  the  aerosol  particles.  Because  aerosol  particles  are  considerably  larger  than 
air  molecules,  aerosol  backscatter  is  typically  much  less  sensitive  to  wavelength.  The 
aerosol  scattering  wavelength  dependence  is  modeled  as  X'x,  where  x  is  the  Angstrom 
coefficient  (x  typically  varies  between  1  and  2). 


Table  1.  Particles  responsible  for  atmospheric  scattering. 


Type 

Radius  (pm) 

Concentration  (cm'3) 

Air  molecules 

io-4 

1019 

Aitken  nucleus1 

10'3  -10'2 

o 

1 

<N 

o 

Haze  particles2 

10'2-1 

•  10-103 

Fog  droplets 

1-10 

10-100 

Cloud  droplets 

MO 

10-300 

Raindrops 

102-104 

10'5  -10‘2 

‘Salt,  dust,  pollen,  etc. 
2Small  fog  droplet. 


Both  molecular  and  aerosol  scattering  contribute  to  the  atmospheric  extinction  (km*1) 
and  backscatter  ( knT'sr'1)  coefficients 

a(z)  =  crR(z)  +  oa(z),  (1) 

P(z)  =  Pr(z)  +  pa(z), 

where  crR(z)  ,pR(z)  and  cra(z)  ,pa(z)  are  the  Rayleigh  (molecular)  and  aerosol 
extinction  /  backscattering  coefficients,  respectively. 

The  aerosol  backscatter  coefficient  profile  gives  a  direct  indication  of  aerosol  loading 
and  qualitatively  reflects  the  shape  of  the  mass-concentration  profile,  whereas  the 
surface  -area  -concentration  profile  is  correlated  with  the  extinction  profile. 

The  backscatter  ratio  ( Rs )  is  commonly  used  to  characterize  the  relative  contributions 
of  molecular  and  aerosol  scattering.  Rs  is  defined  as  the  total  backscatter  coefficient 
divided  by  the  molecular  backscatter  coefficient: 

R  (Z)  =  =  1  +  A  .  (2) 

Pr  Pr 

The  backscatter  ratio  is  derived  from  the  lidar  profile  of  backscattered  energy. 

The  molecular  backscatter  coefficient  profile  can  be  computed  from  theory  using 
model  atmosphere  values  for  the  temperature  T(z)  (degrees  Kelvin)  and  pressure  p(z) 
(millibars): 

/ 

fi.M-  3.6xio-,,.r,#11’.^.  (3)  : 

T(z ) 

The  backscattering  is  optimal  if  the  wavelength  is  of  the  order  of  magnitude  of  the 
mean  particle  radius.  Thus  at  high  altitudes  the  355  nm  return  signal  is  much  larger 
than  the  532  and  1064  nm  signals,  even  though  the  355  nm  laser  pulse  energy  is 
smaller.  Conversely,  the  355  nm  signal  experiences  more  attenuation,  and  the  longer 
wavelengths  penetrate  better  to  the  lowest  layers  of  the  atmosphere.  The  wavelength 
dependence  of  the  particle  backscattering  generally  decreases  with  increasing  particle 
size.  In  [1]  shown  that  the  ratio  of  the  backscattering  at  532  nm  and  355  nm  can  be 


used  to  discriminate  between  aerosol  particles  with  different  size  in  the  upper 
troposphere.  If  a  layer  is  composed  of  aerosols,  then  the  ratio  P532  /  P355  is 
significantly  less  than  unity.  If  the  layer  is  a  thin  cloud  composed  of  ice  crystals, 
which  are  much  larger  than  aerosol  particles,  then  P532  /  P355  is  close  to  unity. 

Table  2  shows  the  backscattering  ratios  for  different  laser  wavelengths  that  were 
calculated  for  a  visibility  of  10  km  using  a  continental  aerosol  distribution. 


Table  2.  Ratio  of  aerosol-to-molecular  scattering. 


Wavelength  [pm] 

Ratio 

0.53 

2.8 

1.06 

23 

1.54 

70 

10.6 

>2000 

Another  valuable  aerosol  quantity  is  the  lidar  ratio  (Rl),  defined  as  the  ratio  of 
aerosol-extinction  coefficient  to  aerosol-backscatter  coefficient  at  a  given  wavelength. 
This  quantity  can  be  used  to  locate  regions  of  atmospheric  aerosol  layering  and  to 
provide  information  on  the  transmission  and  reflection  properties  of  the  atmosphere.  It 
is  a  relative  measure  of  aerosol  size  distribution  with  respect  to  height,  with  larger 
particles  generally  being  represented  by  smaller  lidar  ratios.  For  instance,  a  change  of 
the  lidar  ratio  from  18  to  36  sr  corresponds  to  a  change  of  the  effective  radius  from 
0.35  to  0.18  pm. 

The  actual  values  of  Rl  for  tropospheric  aerosol  may  range  from  28  to  48  sr  (typical 
of  rural  atmosphere).  In  [2]  were  suggested  values  for  the  lidar  ratio  of  30  sr  in  the  , 
troposphere  and  50  sr  in  the  stratosphere.  In  [3]  is  given  a  lidar  ratio  of  about  60  sr  for 
small  stratospheric  liquid  droplets  with  an  effective  radius  of  about  0. 1  pm  and  of 
about  10  sr  for  cirrus  ice  crystals  (10  to  1000  pm  in  size). 

Radiative  and  microphysical  properties  of  aerosols  are  variable  in  time  and  space. 

The  planetary  boundary  layer  (PBL,  2-3  km)  contains  most  of  the  aerosol  and  water 
vapor  in  the  atmosphere.  Aerosols  in  the  PBL  generally  have  the  accumulation 
particle  mode  around  0.1 -0.2  pm  in  radius  and  the  large  particle  mode  around  several 
micrometers  for  volume  size  distribution.  The  large  mode  consists  mainly  of 
soil-derived  particles  or  sea  salt  particles  whereas  the  accumulation  mode  consists  of 


anthropogenic  particles  such  as  sulfates  or  nitrates.  Maximal  values  of  the  individual 
modes  change  with  season  or  with  different  atmospheric  conditions. 

Aerosols  in  the  free  troposphere  (4-12  km)  are  from  two  main  sources:  aftermaths  of 
major  volcanic  eruptions  and  desertic  dust  episodes.  Studies  of  Saharan  dust  optical 
properties  show  that  the  wavelength  dependence  of  the  aerosol  extinction  coefficient 
is  almost  zero  in  the  visible  and  the  near  UV  owing  to  an  enhancement  of  the  number 
of  large  diameter  particles  (>2  pm). 

The  base  of  the  stratospheric  aerosol  layer  is  usually  at  the  tropopause  (approximately 
16  km).  The  main  structure  of  the  stratospheric  aerosol  profile,  the  Junge  layer,  is 
located  a  few  km  above  the  tropopause.  Particles  greater  than  a  micrometer  in  size  are 
more  effectively  removed  from  the  stratosphere  because  of  the  Earth's  gravitation. 
Above  the  stratospheric  aerosol  layer,  the  laser  backscatter  becomes  purely  Rayleigh 
and  can  be  related  directly  to  molecular  density. 

Aerosol  profiles  are  summarized  in  Table  3  [4]. 

The  aerosol  concentration  is  expressed  in  terms  of  extinction  coefficients  for  a 
reference  wavelength  of  0.5  pm.  Depending  on  the  aerosol  type  three  to  four  different 
layers  are  considered.  The  first  and  second  layers  describe  the  distribution  of  aerosol 
in  the  troposphere.  The  third  and  the  fourth  layers  indicate  the  vertical  state  of  the 
tropospheric  background  aerosol  and  the  stratospheric  aerosol. 

Maritime-mineral  aerosols  are  considered  differently,  as  the  layer  up  to  4.5  km  is 
divided  into  two  layers.  The  first  layer  consists  of  pure  maritime  aerosols,  while  the 
second  layer  is  populated  by  freshly  produced  dust.  The  extinction  coefficient  values 
in  Table  3  refer  to  dry  particles. 

The  aerosol  size  distribution  for  all  types  is  given  in  [4,5;6].  and  in  a  previous  report  . 
(December  1998). 


Table  3  .  Vertical  distribution  of  different  atmospheric  aerosol  types. 

<yeo  and  oe0  (h)  give  the  extinction  coefficient  in  km'1  for  10  m  above  ground  and  at  a 
height  h  respectively.  The  quantity  Z  means  the  scale  height  for  an  exponential  profile 


type. 


Aerosol  type 

1. layer  2. layer 

3. layer 

4. layer 

(trop.  backgr.) 

(stratosph.) 

clean  continental 
(rural) 

0-2  km 

CTe(h)=CTeo  =0.126 

2-12  km 
0.02-0.005 

12-35  km 
0.01-2E-4 

average  continental 

0-2  km 

ae(h)=<Teo  =0.22 

2-12  km 
0.02-0.005 

12-35  km 
0.01-2E-4 

desert 

(wintertime) 

0-3  km 

<7e(h)=<7eo  =0-303 

3-12  km 
0.02-0.005 

12-35  km 
0.01-2E-4 

desert 

(summertime) 

0-3  km 

CTe(h)=aeo  =0.735 

3-12  km 
0.025-0.005 

12-35  km 
0.012-2. 5E-4 

urban 

0-2  km 

ae(h)=CTeo  =0.436 

2-12  km 
0.02-0.008 

12-35  km 
0.01-2.2E-4 

clean  maritime 

0-3  km 

ae(h)=  CTeo  e'h/z 

CTeo  =0.123, 

Z=1  km 

3-12  km 
0.02-0.002 

12-35  km 
0.005-2E-4 

maritime  mineral 

0-3  km  3-4.5  km 

CTe(h)=CTeo  e'h/Z  ,  CTeo  =0.342 
CTeo  =0.125, 

Z=1  km 

4.5 -12  km 
0.03-0.005 

12-35  km 
0.008-2E-4 

9 

maritime  polluted 

0-3  km 

CTe(h)=CTe<re'h/z 

3-12  km 
0.012-0.005 

12-35  km 
0.008-2E-4 

CTeo  =0.147, 

Z=lkm 

Arctic  aerosol 

0-3.5  km 

CTe(h)=CTeo  e'h/Z  , 

CTeo  =0.317,  Z=1.4  km 

3.5-12  km 
0.015-0.005 

12-35  km 
0.008-2E-4 

Antarctic  aerosol  0-10  km  10-12  km  12-35  km 

<7e(h)=  Geo  e'h/z ,  0.002  2.2E-4 

CTeo  =0.018,  Z=3.5  km 

Antarctic  aerosol 


10-12  km 

0.002 


12-35  km 
2.2E-4 


2.  Atmospheric  transmission  and  optical  depth. 


The  laser  irradiance  in  the  target  plane  located  at  range  z  from  the  transmitter  is  given 
by 


nW\z)  ’ 


(4) 


w(z)=fff7¥~7,  (5) 

where  W(z)  represents  the  radius  of  the  laser  pulse  in  the  target  plane  at  the  instant  of 
interest,  Wo  is  the  laser  output  aperture  radius,  0  is  half  the  laser  beam's  divergence 
angle.  Pi  is  the  total  output  power  of  the  laser,  and  T(X,z)  is  the  atmospheric 
transmission  at  laser  wavelength  to  range  z .  The  distribution  of  laser  irradiance  across 
the  target  plane  is  assumed  to  be  uniform  over  the  area  of  illumination. 

The  transmission  or  attenuation  factor  T(X,z)  follows  the  exponential  law  of 
attenuation,  the  Beer-Lambert  law.  It  is  related  to  the  path  optical  depth  or  integrated 
extinction  by 

T(X,z)  =  exp(-r(A,  z))  (6) 

or 

z 

T(X,z)  =  Qxy(-\cr(X,z)dz),  0) 

0 

where  r(z)  is  the  atmospheric  optical  depth  through  height  z  and  a(z)  is  the 
atmospheric  unit  volume  extinction  coefficient  (km"1)  at  any  intervening  height  z 
between  0  and  z.  The  extinction  coefficient  includes  the  loss  effects  of  both  scattering 
and  absorbtion. 


In  Figs.2-7  the  horizontal  profiles  of  optical  depth,  atmospheric  transmission  and  laser 
irradiance  in  the  target  plane  for  different  aerosol  types  at  different  elevation  heights 
are  presented. 


Fig.2.  Horizontal  profiles  of  optical  depth  at  different  elevation  heights.  Desert 
aerosol  type. 


Fig.3  Horizontal  profiles  of  atmospheric  transmission  at  different  elevation  heights 
Desert  aerosol  type. 


Fig.  4.  Horizontal  profiles  of  optical  depth  at'different  elevation  heights.  Average 
continental  aerosol  type. 


Fig.  5.  Horizontal  profiles  of  atmospheric  transmission  at  different  elevation  heights 
Average  continental  aerosol  type. 


Fig.  6.  Horizontal  profiles  of  laser  irradiance  (Watt  /  m2 )  in  the  target  plane  at 
different  elevation  heights  for  transmitter  energy  of  1  J  (7  ns  pulse).  Average 
continental  aerosol  type. 


Fig.7.  Horizontal  profiles  of  laser  irradiance  (Watt  /  m2 )  in  the  target  plane  at 
different  elevation  heights  for  transmitter  energy  of  1  J  (7  ns  pulse).  Desert  aerosol 
type. 


3.  Laser  beam  spreading  in  the  atmosphere. 


The  effect  of  the  atmosphere  on  the  intensity  of  the  laser  beam  in  the  target  plane  is  to 
spread  the  laser  beam  out  because  of  turbulent  and  aerosol  particle  scattering  and  thus 
to  reduce  the  total  energy  flux  arriving  at  the  target  plane. 

The  average  irradiance  <I>  (watt  /  cm2)  at  a  distance  p  from  the  optical  axis  of  the 
transmitter  in  a  target  plane  located  at  range  z  is  given  by  [8] 

<I(p,z)>*ex p[-<ra  •  z]  ■  exp[- p2  / W2 ]+  ~~p ' exp[- P2  !Ws  \  >  (8) 

7VWt  7tW  s 

cra-z»  1 , 
where 


W2  =  W2  + 

rr  *  vac 


\kPtj 

W2=w2+w2, 

W-2  _  0*0  '  2  •  Qrms 

na  ~  -  ’ 


(9) 

(10) 
(ID 


w;ac  =  w2 


i  + 


kW, 


0  J 


pt  =  (o.5  •  k 2  ■  C2  ■  z ) 3/5  for  spherical  waves,  ca  is  the  aerosol  scattering  coefficient, 
k  =  2n  /  X . 

02ms  =<  <9 2  >  is  the  mean  squared  scattering  angle  averaged  over  all  the  aerosols  [6]. 

Wt  and  Wa  give  the  contribution  to  the  beam  width  due  to  the  turbulence  and  aerosol 
scattering  respectively.  Wvac  is  the  beam  spot  size  in  the  vacuum. 

The  aerosol  absorption  and  molecular  scattering  effects  are  assumed  negligibly  small 
at  selected  wavelengths  (>0.3  pm). 

For  a  laser  beam  to  detect  a  species  and  targets  at  a  long  distance,  the  beam  must  not 
be  appreciably  attenuated  by  the  intervening  atmosphere.  The  output  wavelength  of 
the  laser  must  therefore  lie  in  a  spectral  transmission  window  of  the  atmosphere.  The 


most  useful  transparent  spectral  ranges  are  the  visible  (0.4-0.7  pm),  near-infrared 
(0.7-1. 5  pm),  and  3  to  5  pm  and  9  to  13  pm  regions. 

At  the  Nd:YAG  wavelengths  (1064  and  532  nm)  gas  absorption  becomes  negligible 
as  these  wavelengths  are  located  in  a  propagation  window  where  no  local  line 
absorption  structure  or  continuum  absorption  bands  exist. 

In  the  first  term  of  Eq.  8,  the  aerosol  scattering  extincts  the  intensity  while  turbulence 
spreads  out  the  beam  from  its  vacuum  value  according  to  Eq.  9.  In  the  second  term, 
scattering  from  aerosols  further  expands  the  beam. 

Ref.  [9]  gives  the  average  effective  beam  spot  size  (beam  radius)  for  a  beam  with  a 
spot  size  in  free  space,  including  beam  spread  and  beam  wander  for  the  propagation 
path  through  turbulence  without  the  aerosol  contribution.  It  is  derived  from  the 
average  on-axis  intensity  and  is  given  by 

W,  =fP„$  +  6.72*,,‘  •  K-ml”  \(H -hJIax&Y'V ’  (12) 


H 

where  m0  =  jc2(h)dh ,  ho  is  the  laser  source  altitude  ,  H  is  the  target  altitude,  £is  the 

fy) 

zenith  angle,  A  =  2z/kWvac2 ,  where  Wvac  is  the  diffraction-limited  beam  radius  at  the 
target  plane,  and  C2n  is  the  turbulence  structure  constant  (m'2/3). 

Also  presented  there  is  the  average  beam  spot  size  that  includes  beam  spread  but  not 
beam  wander,  i.e.,  the  physical  radius  of  the  effective  propagating  beam, 


W  = 

"  n 


f  1  ] 

1/3 

W,  -1.96 Wl  ■ 

•  k  ■  z  •  A  •  m0 

K2o}o) 

(13) 


where  cd0  is  the  transmitted  beam  radius. 

In  Figs.  8-21  expected  laser  beam  spot  size  versus  distance  in  vacuum  with  the 
turbulence  and  particle  scattering  contributions  for  different  aerosol  types  and 
turbulence  strengths  are  presented.  The  turbulence  strength  vertical  profiles  are  given 
in  a  previous  report  (December  1998). 

In  these  calculations  the  laser  wavelength  is  assumed  to  be  532  nm.  In  general  it 
was  found  that  the  laser  beam  widening  caused  by  aerosols  is  very  significant, 


even  more  so  and  sometimes  considerably  more  so,  than  that  caused  by 
turbulence.  However,  the  aerosol  widening  is  very  wavelength  dependent.  When 
effective  particle  radius  is  larger  than  wavelength,  laser  beam  widening  caused 
by  aerosols  diminishes. 


Fig.  8.  Laser  beam  spot  size  (m)  as  a  function  of  horizontal  distance  at  2  km  altitude. 
Average  continental  aerosol  type.  Cn2=l.le-16  m‘2/3  (daytime  SLC  model). 


103  104  10s  z,m 

Fig.  9.  Laser  beam  spot  size  (m)  as  a  function  of  horizontal  distance  at  10  km  altitude. 
Average  continental  aerosol  type.  Cn2=2e-18  m'2/3  (daytime  SLC  model). 
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Fig.  12.  Laser  beam  spot  size  (m)  as  a  function  of  horizontal  distance  at  10  km 
altitude.  Desert  aerosol  type.  C„2=2e-18  m'2/3  (daytime  SLC  model). 
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Fig.  13.  Laser  beam  spot  size  (m)  as  a  function  of  horizontal  distance  at  20  km 
altitude.Desert  aerosol  type.  Cn2=1.41e-18  m’2/3  (daytime  SLC  model). 


Fig.  14.  Ratio  of  beam  spot  size  with  turbulence  widening  compared  to  that  in  vacuum 
as  function  of  horizontal  distance  at  10  km  altitude.  Desert  aerosol  type.  Cn2  =  2e-18 
m'2/3  (daytime  SLC  model). 
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Fig.  15.  Ratio  of  beam  spot  size  with  aerosol  widening  compared  to  turbulence 
widening  as  function  of  horizontal  distance  at  10  km  altitude.  Desert  aerosol  type.  Cn2 
-  2e-18  m'2/3  (daytime  SLC  model). 


Fig.  16.  Ratio  of  beam  spot  size  with  turbulence  widening  compared  to  that  in 
vacuum  versus  horizontal  distance  at  10  kmaltitude.  Average  continental  aerosol  type. 
Cn2  =  2e-18  m'2/3  (daytime  SLC  model). 


Fig.  17.  Ratio  of  beam  spot  size  with  aerosol  widening  compared  to  turbulence 
widening  as  function  of  horizontal  distance  at  10  km  altitude.  Average  continental 
aerosol  type.  C„2  =  2e-18  m‘2/3  (daytime  SLC  model). 


Fig.  18.  Ratio  of  beam  spot  size  with  turbulence  widening  compared  t  that  in  vacuum 
as  function  of  horizontal  distance  at  20  km  altitude.  Desert  aerosol  type.  Cn2  = 

1 .4 1  e- 1 8  m'2/3  (daytime  SLC  model). 
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Fig.  19.  Ratio  of  beam  spot  size  with  aerosol  widening  compared  to  turbulence 
widening  as  function  of  horizontal  distance  at  20  km  altitude.  Desert  aerosol  type. 
C„2  =  1.41  e- 18  m'2/3  (daytime  SLC  model). 


Fig.  20.  Ratio  of  beam  spot  size  with  turbulence  widening  compared  to  that  in  vacuum 


versus  horizontal  distance  at  20  km  altitude.  Average  continental  aerosol  type. 


Fig.  21.  Ratio  of  beam  spot  size  with  aerosol  widening  compared  to  turbulence 
widening  as  function  of  horizontal  distance  at  20  km  altitude.  Average  continental 
aerosol  type.  Cn2  =  1.41  e- 18  m'2/3  (daytime  SLC  model). 


4.  Lidar  system  discription. 


The  lidar  we  are  developing  to  measure  vertical  profiles  of  aerosol  size  distribution 
and  turbulence  strength  is  almost  ready.  All  the  main  components  have  been  derived 
upon  and  ordered.  It  should  be  operational  by  the  end  of  the  year. 

Fig.22  shows  the  schematic  of  the  lidar  setup,  composed  of  a  laser,  receiving 
telescope,  signal  processing  unit,  and  data  processor.  Table  4  gives  the  specifications 
of  the  various  subsystems. 

The  lidar  system  employs  aNd:YAG  laser  operating  on  the  fundamental  (1064  nm), 
doubled  (532  nm),  and  tripled  (355  nm)  wavelengths.  The  receiving  telescope  is  a 
Classical  Cassegrain  design  with  a  diameter  of  12.5  inch.  The  optical  subsystem  in  the 
telescope  uses  dichroic  beam  splitters  to  direct  the  three  wavelengths  to  different 
detectors.  The  532-  and  355-nm  channels  are  detected  with  photomultiplier  tubes 
(PMT)  operating  in  analog  mode,  and  the  1064  nm  channel  is  detected  with  a  silicon 
avalanche  photodiode  (APD).  Detected  returns  at  all  three  wavelengths  are  to  be  used 
to  calculate  aerosol  size  distribution  at  each  elevation  range.  Vertical  resolution  is 
expected  to  be  15  m.  In  addition,  a  cooled  Intensified  CCD  camera  (ICCD,  Andor 
Technology)  used  in  a  novel  mode  as  the  receiver  for  extracting  the  atmospheric 
turbulence  profiles.  The  ICCD  has  a  time  gating  and  delay  that  allowing  to  obtain  a 
series  of  time-  and  altitude-resolved  signals. 

Radiation  from  the  Nd:YAG  laser  is  directed  into  the  atmosphere  via  a  right-angle 
prism. 

The  transmitting  optics  and  the  telescope  are  separated  by  approximately  1 .2  m  with 
respect  to  their  optical  axes  to  form  a  noncoaxial  system.  Hence  the  backscattered 
echoes  are  not  seen  by  the  telescope  until  the  transmitted  pulse  crosses  into  and  starts 
overlapping  with  the  telescope  field  of  view.  This  reduces  the  signal-induced  noise  of 
the  PMT  (induced  by  the  strong-intensity  signal  from  low-altitude  regions,  which  can 
get  superimposed  with  the  weak  signals  from  the  high  altitudes). 

The  overlap  range  may  be  extended  for  long-range  probing  applications  (stratospheric 
sounding)  to  avoid  near  range  signal  saturation  of  the  photodetector.  Alternatively, 
mechanical  shutters  may  be  used  to  block  light  from  reaching  the  detector  until  after 
some  delay.  Narrow-band  interference  filters  (~1  to  2  nm  bandwidths)  are  also 
included  in  the  receiver  system  to  reduce  background  radiation. 


The  transient  signal  response  to  each  transmitted  lidar  pulse  is  sampled  and  digitized. 
Sampling  rate  is  on  the  order  of  one  sample  per  100  ns  (15  m  resolution).  Each 
sample  is  digitized  with  14-bit  resolution.  To  improve  the  signal-to-noise  ratio,  the 
digitized  signals  for  each  transmitted  lidar  pulse  are  averaged  over  tens  to  hundreds  of 
pulses. 


Fig.  22  Block  diagram  of  the  lidar  system. 


Transmitter 

Laser 

Output  energy/pulse 


Repetition  rate 
Pulse  duration 
Beam  divergence 


Receiving  optics 


Telescope 
Field  of  view 


Nd:YAG,  Continuum  model  Surelite  I- 10  with  2  harmonics 
450  mJ  at  1064  nm 
200  mJ  at  532  nm 
60  mJ  at  355  nm 
10  Hz 
7  ns 

0.6  mrad 

12.5”  Classical  Cassegrain  (Davro  Optical  Systems,  Inc.) 
2mrad 


Filter  linewidth 


lnm 


Filter  transmittance  45%  (1064nm  and  532nm) 

25%  (355nm) 

Detectors 

Intensified  CCD  camera  Andor  Technology  Inc. 

Photomultiplier  tube  Hamamatsu  R4632 

Quantum  efficiency  25%  at  355nm 


18%  at  532nm 


Current  gain 
Dark  current 
Excess  noise  factor 
Avalanche  photodiode 
APD  quantum  efficiency 
APD  responsivity 
APD  excess  noise  factor 
APD  dark  current 


Signal  processor 
Type 

Minimum  sample  rate 


Accuracy 

Elevation  resolution 


3.5xl06 

InA 

1.2 

EG&G  C30954E 
30%atl064nm 
36A/W 
4 

50nA 

Digital  processing 
100  nsec 

14  bits 

15  m 


Table  4.  Lidar  System  Specifications 


The  lidar  response  may  be  interpreted  in  terms  of  the  lidar  range  equation. 

The  instantaneous  received  power,  P(z),  due  to  backscattering  from  height  z, 
assuming  only  single  scattering  and  a  vertically  pointing  lidar,  may  be  expressed  in 
the  form: 


P(z)  =  P0-G(z) -m-T'iz)  ,  (13) 

2  z 

where 

Po  =  transmitted  peak  power, 

G(z)  =  geometry  factor  for  overlap  of  transmitter/receiver  beam  paths, 

ctp  =  transmitted  pulse  length  (tp  is  the  pulse  duration,  c  is  the  speed  of  light), 

Ar/ z2  =  solid  angle  (sr)  subtended  by  receiver  aperture  Ar  at  range  z, 

P(z)  =  unit  volume  backscattering  coefficient  (rn'sr'1)  of  the  atmosphere  at  range  z, 
^(z)  =  atmospheric  transmission  to  range  z  and  back  to  the  receiver. 

The  geometric  form  factor  can  be  evaluated  fairly  easily  for  a  biaxial  lidar  (Fig.23)  if 
the  objective  lens  (or  mirror)  of  the  telescope  represents  the  limiting  aperture  of  the 
receiver  optics,  any  obstruction  is  neglected,  and  an  uniform  laser  distribution  over 
the  area  of  illumination  is  assumed. 


Fig.  23.  Geometry  of  biaxial  lidar. 


Under  the  conditions  of  Fig.23  one  can  think  of  G(z)  as  a  simple  overlap  factor,  and 
one  can  write 


(z),r(z);rf(z)} 
1  }  nW\z) 


(14) 


where  A  represents  the  area  overlap  function,  rx(z)  is  the  radius  of  the  receiver  optics 
field  of  view  in  the  target  plane  at  range  z,  and  W(z)  is  the  radius  of  the  circular 
region  of  laser  illumination  or  spot  at  range  z.  In  formulating  these  equations,  the  field 
of  view  of  the  receiver  optics  in  the  target  plane  is  assumed  a  circle  of  radius 


rT(z)  =  r0  +(j>-z. 


(15) 


where  ro  is  the  effective  radius  of  the  telescope  lens  (or  mirror),  and  <|>  is  the  receiver 
optics  half  opening  angle. 

As  for  the  laser  beam  radius,  assuming  the  output  of  the  laser  is  predominantly  in  the 
TEMoo  mode,  one  can  write 

w(  2)=ylwl7er^7,  (i6) 

where  Wo  is  the  laser  output  aperture  radius,  and  0  is  the  laser's  half  divergence  angle. 
The  separation  of  the  telescope  and  laser  axes  in  the  target  plane  is 

d  =  d0-z-S,  (17) 


where  do  is  the  separation  of  the  axes  at  the  lidar,  and  5  is  the  inclination  angle 
between  the  laser  and  the  telescope  axes. 

Three  situations  are  possible: 

1 .  The  separation  of  the  axes  is  small  enough  that  either  the  area  of  laser  illumination 
lies  totally  within  the  receiver  optics  field-of-view  and  correspond  to  an  overlap 
factor  equal  to  unity.  This  amounts  to  saying  that  if  d  <  | rT  -  W |  than  the  overlap 


factor  may  be  computed  as  the  ratio  of  areas.  That  is  G(z)  = 


*7y20) 

7tW2(z) 


(l)2 


2.  The  separation  of  the  axes  is  too  large  for  there  to  be  any  overlap  between  the 
receiver  optics  field-of-view  and  the  area  of  laser  illumination,  that  is,  A=0  if 
d>  rr  +W. 

3.  The  separation  of  the  axes  lies  between  these  extremes: 

\rT-W\<  d<rT  +W. 

Under  these  circumstances  the  area  overlap  function  is 

A{rT,W,d}  =  W2  -y/w  +r 2  -y/r  -rT  •^•sin(^'r),  (18) 


where 


Vw 


=  arccosl 


</2  +  lf2  -a>2 
2  Wd 


Vr 


=  arccosl 


d2 


+  rr2  -  IT2 
2  rTd 


(19) 


Clearly,  the  overlap  factor  G(z)  will  depend  on  the  range  z  and  on  the  values  of  the 
angular  parameters  5,  9,  and  (f>. 

The  feature  of  a  biaxial  arrangement  of  the  lidar  system  cam  be  useful  in  eliminating 
near-field  scattering  -  an  effect  that  can  saturate  the  photodetector  when  the  lidar  is 
designed  for  long-range  measurements.  Thus,  for  the  laser  and  telescope  axes  being 
parallel  (i.e.,  8=  0),  the  variations  of  G(z)  with  z  for  three  values  of  do  (0.8  m,  1  m, 
and  1.2  m)  are  presented  in  Fig.  24. 

The  behavior  of  overlap  factor  G(z)  with  changing  half  the  opening  angle  of  the 
telescope,^,  are  presented  in  Fig.  25,  where  ^varies  between  1  and  3  mrad.  The  laser 
divergence  half  angle  9=  0.3  mrad. 

The  behavior  of  G(z)  with  an  inclination  between  the  telescope  and  laser  axes  is 
presented  in  Fig.  26 


Fig.24.  Variation  of  overlap  factor  G(z)  with  range  z  for  different  separation  of  axes  at 
lidar  do .  The  receiver  half  opening  angle  ^=lmrad  and  3=0. 


Fig.25.  Variation  of  overlap  factor  G(z)  with  z  for  the  different  receiver  field-of-view 
do=  1 .2  m  and  3=0. 


Fig.  26.  Variation  of  G(z)  with  range  z  for  different  inclination  between  telescope  and 
laser,  do  =  1.2  m  and  d>=  lmrad. 


These  results  clearly  indicate  that  interpretation  of  short-range  lidar  measurements 
must  take  info  account  the  geometrical  factors  involved. 

For  measurements  at  shorter  ranges  the  axes  of  the  transmitter  and  receiver 
field-of-view  can  be  made  to  intersect  at  some  finite  height  down  to  ~600  m.  In  this 
case  the  neutral  density  filter  must  be  used  in  front  of  the  detector  to  reduce  the  strong 
near-field  radiation.  Also,  the  detector  gain  can  be  decreased. 

For  the  simulation  of  the  backscattering  lidar  signals,  two  extinction  profiles  are 
selected  that  are  considered  possible  cases  for  aerosol-polluted  cloud-free 
atmospheres  [10].  The  profiles  describe  a  "clear"  and  "hazy"  atmosphere 
corresponding  to  visibilities  of  23  and  5  km  respectively  at  ground  level.  The  vertical 
profiles  of  the  aerosol  extinction  coefficient  at  532  nm  are  shown  in  Fig.27. 

Fig  28  shows  the  backscattered  single-shot  signal  for  532  nm  wavelength.  The 
sharply  rising  initial  portion  of  the  signal  is  due  to  the  geometrical  overlap  of  the 
transmitted  beam  and  the  receiver  field  of  view.  Following  system  convergence  the 
signal  decreases  as  1/  R2  times  the  Beer's  law” factor. 

The  lidar  system  performance  can  be  estimated  by  deriving  the  signal-to-noise  ratio  as 
a  function  of  altitude. 

A  simple  signal-to-noise  ratio  equation  for  a  lidar  using  incoherent  detection  is 


S  K 

N  JN,+Nt+N,  ’ 


(20) 


where  Ns  =  P(z)  •  tj  •  At  ■  JU  he  is  the  number  of  signal  photons  detected  in  sample 

time  At,  Nb  is  the  number  of  detected  photons  resulting  from  the  background 
illumination  in  time  At,  Ne  represents  the  noise  current  of  the  detector,  and  q  is  the’ 
detector  quantum  efficiency.  For  avalanche  photodiode  detector  Ne  is  often  the 
dominant  noise  source.  Since  a  PMT  detector  offers  high,  nearly  noise-free  gain,  Ne  is 
usually  negligibly  small  in  that  case.  During  night  operation,  Nb  is  also  often 
negligibly  small.  In  this  case  the  dominant  noise  source  is  the  statistical  noise  of  the 

signal  itself  and  SNR  oc  .  Fig.  29  shows  the  SNR  versus  altitude. 

The  goal  of  lidar  inversion  is  to  retrieve  the  optical  atmospheric  parameters  a  and  (5 
from  the  return  signal  P(z).  A  number  of  algorithms  were  presented'  in  the  previous 
report  (September  1998). 


Fig.  28.  Backscattered  single-shot  lidar  signal  (Watt)  as  a  function  of  range 
A.-532  nm,  "Clear"  atmosphere. 


Fig.  29.  Signal-to-noise  ratio  as  function  of  altitude.  X=532  nm. 


5.  Conclusion. 


Laser  beam  broadening  for  elevations  studied  here  (2-20  km)  results 
significantly  from  aerosols,  usually  even  more  so  than  by  turbulence,  as  shown  in 
Section  3  here. 
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